Localization and absence of Breit-Wigner form for Cauchy random band matrices 
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We analytically calculate the local density of states for Cauchy random band matrices with 
strongly fluctuating diagonal elements. The Breit-Wigner form for ordinary band matrices is re- 
placed by a Levy distribution of index n = 1/2 and the characteristic energy scale a is strongly 
enhanced as compared to the Breit-Wigner width. The unperturbed eigenstates decay according 



to the non-exponential law cx e 



We analytically determine the localization length by a new 



method to derive the supersymmetric non-linear a model for this type of band matrices. 



PACS numbers: 72.15.Rn, 05.40. Fb, 05.60.Gg 

Band matrices with random elements appear in a var- 
ity of physical problems in the context of quantum chaos 
and localization A detailed analytical investiga- 

tion of the properties of such matrices was performed 
by Fyodorov and Mirlin Q using Efetov's supersymme- 
try technique Q. Later, motivated by the localization 
problem of two interacting particles, Shepclyansky intro- 
duced and studied random band matrices superimposed 
with strongly fluctuating diagonal matrix elements 
Subsequent work on this type of matrices (jj showed the 
existence of the Breit-Wigner regime where the eigen- 
states have a very peaked structure inside an eventual 
localization domaine. 

In all these cases the matrix elements were drawn from 
a regular, typically gaussian, distribution with finite vari- 
ance. In the present work, we extend the band matrix 
ensembles of refs. [||j7j to the case of Cauchy distributed 
matrix elements. In an early work, Lloyd J8| already in- 
troduced and studied a model with diagonal Cauchy dis- 
tributed disorder. A detailed study of the level statistics 
of full random matrices H with elements drawn from 



the more general Levy distribution |10| was performed 
by Cizeau and Bouchaud. Band matrix ensembles with 
Cauchy distributed matrix elements were recently argued 
to be relevant and studied numerically in the context 
of the localization problem of two interacting particles 
[|ll],[l2|]. In this work, we present analytical results con- 
cerning the wave function properties for a model similar 
as in ref. jl2| but which differs from that of ref. || by 
two important features: it concerns a banded and not 
full matrix, and its diagonal elements are typically much 
larger than the off-diagonal coupling elements. 

Actually, we think that this case is of particular in- 
terest for a generic disordered fermionic many particle 
problem (for instance recently studied in p3[ ) where non- 
interacting eigenstates are coupled by quasi-random two- 
particle interaction matrix elements. There is an interest- 
ing and subtle regime jl4j where the effective two-particle 
level spacing is smaller than the typical interaction ma- 
trix element while typical eigenstates are nevertheless 
composed of many non-interacting eigenstates. In this 
regime, we can diagonalize the problem for two parti- 



cles by perturbation theory and apply the corresponding 
unitary transformation to the full many-particle Hamil- 
tonian. This provides a new type of random matrix 
with residual three-body interaction matrix elements typ- 
ically given by perturbative expressions with denomina- 
tors, containing independent random non-interacting en- 
ergy differences. Iterating this procedure one also ob- 
tains higher order contributions. Flambaum et al. 
have argued that such denominators give effective distri- 
butions for these matrix elements characterized by long 
tails with the same power law as the Cauchy distribution. 
Even though the real situation is more complicated, this 
argument indeed shows the relevance of Cauchy random 
matrix ensembles. 

The model we consider is a random band matrix of 
dimension N x N and of width 6^1 whose matrix ele- 
ments are of the form H^i = i]k Ski +Uki with Uki = if 
k = I or \k — l\ > b. These matrix elements are statisti- 
cally independent and distributed according to a Cauchy 
distribution p a (x) = n^ 1 a/(a 2 + x 2 ) where the width 
a is given in terms of two parameters W and Uq via 
a = W/ir for r]k and a = Uq for Uu- In this paper, we 
consider the orthogonal symmetry class where H is real 
symmetric. The generalization to the other symmetry 
classes is straightforward. For the following analytical 
calculations, we will furthermore concentrate on the case 
where Uq is sufficiently large to avoid a simple perturba- 
tive situation but still so small that the total density of 
states will essentially be dominated by the unperturbed 
diagonal elements only. 

For ordinary band matrices |Q with a finite variance 
of Ujk, this case corresponds to the Breit-Wigner regime 
where the eigenstates are in addition to an eventual space 
localization also localized in the (unperturbed) energy 
space, i. e. only the sites j with \E — rjj\ < T es- 
sentially contribute to an eigenstate of energy E where 
T w 2ir{2b)(\U jk \ 2 )/W < W is the Breit-Wigner width. 
This can be seen in the local density of states at site j 
which is a Lorentzian in (E — r)j) of width T/2 0. 

To study the Cauchy band matrix ensemble introduced 
above, we therefore start with this quantity, 
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Pj (E) = —lm 

7T 



G%\E) 



(1) 



where G^'(E) = (E±iO — H)~ 1 and the average (• • -)j is 
taken with respect to all random matrix elements except 
the diagonal element r/j at the site j under considera- 
tion. As in the original Lloyd model ||] the average over 
the other diagonal elements rjk, k j can be exactly 
performed by replacing in the definition of the Green 
function these elements according to: r/k + i — > i W/n. 
Using a general algebraic identity for the block inverse of 
a matrix, we obtain 



Pj (E) = --Im 

7T 
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E + iQ-r)j + iTj(U)/2 



Tj(U) = 2i 



(2) 
(3) 



where (•••)[/ denotes the average over Ujk only, and 
<§(+) = {E+iW/n-U)- 1 with [/being the (N-l)x(N- 
l)-matrix obtained from U after elimination of the j-th 
row and j-th column. Eq. (^) reproduces the Lorentzian 
Breit-Wigner form of width ~Re\Tj(U)]/2 for a given real- 
ization of the coupling matrix U. For ordinary band ma- 
trices this width is a self-averaging quantity and the fur- 
ther [/-average does not modify the Breit-Wigner form. 
However, for the case of Cauchy band matrices Tj(U) 
is strongly fluctuating and the [/-average will consider- 
ably modify the Breit-Wigner form. Now, we consider 
the limit \E\ <C W and we neglect U in the definition of 
G*( + ) which is possible for the regime we consider. Using, 
(?(+) w (-in/W)l N , we obtain 
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with the function cf>(z) defined by 



4(z) 



du ■ 



1 + u 2 



~ zu ^ 2 , z > . (5) 



In the limit b ^> 1, the integral (Q) is dominated 
by the behavior at small t. Using the approximation 

hi(/)(z) w —\/^z + 0(z), we obtain 



P J {E) = -L 1/ J E Vj 



a \ a 
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Ifi U 2 h 2 

LbU ° b (6) 



W = ± J dt e«-W . 



(7) 



Here L^s) represents a Levy distribution of index /i JTcj ] 
with the behavior L^(s) cx s~ 1_A1 for \s\ ^> 1. The ex- 
pressions (|^) and ([?]) provide the first important results 
of this work. They show that the local density of states is 



still a peaked function of {E — rjj) but with two important 
modifications. First, the Lorentzian Breit-Wigner form 
is replaced by the Levy distribution L1/2, and second, the 
characteristic energy scale behaves as a ~ (4/7r) bT ^> T 
where T is the Breit-Wigner width for ordinary band ma- 
trices (with a finite variance (U 2 k ) = Uq). We mention 
that the approximation to neglect the contribution of U 
in is valid for a < W, i. e. U b < W. 

The first modification has an important implication for 
the time-evolution of a state \ip(t)> initially localized at 
one site, \ip(0) >= \j >■ Then, the average of the ampli- 
tude <ip(t)\j > obeys a non-exponential decay law |ll| 
of the form cx e~ v/ "' instead of cx e~ Tt l 2 . This type of 
behavior was for instance recently found in the context 
of many-body effects in cold Rydberg atoms WA - We 
conclude that it is a very general feature due to the long 
tail distribution of residual interaction matrix elements 
in many-body problems |l5f| . 

The second modification concerning the large enhance- 
ment of the characteristic energy scale can be qualita- 
tively understood by the fact that for a given realization 
of U the sum in eq. (||) is essentially determined by the 
typical maximal value 2b Uq of the 2b different matrix 
elements Ujk- 

It is pretty obvious to generalize eqs. (||) and (]?]) to 

the case where the matrix elements Ujk are drawn from 

the more general Levy distribution Uq 1 L^Ujk/Uo). In 

this case, we have to replace in (||) Lin by L^/ 2 an d 

the characteristic energy scale is enhanced according to: 
a ^&2/M-i r . 

We now turn to a field theoretical formulation of the 
Cauchy band matrix model in terms of a supersymmet- 
ric non- linear a model [Q. This formulation provides 
the access to the correlators of an advanced Green func- 
tion G^ and a retarded Green function G^ which are 
for example useful to study the transport properties, the 
level statistics or the appearance of dynamical localiza- 
tion. For this, typically, the ensemble average of a gen- 
erating functional of the form 

F(J) = Sdet" 1/2 (£+(! + ie)A -H(g>l 8 + J) (8) 

is considered. For the details concerning the super- 
symmetry method and the notations, we refer to the 
standard literature Q^]. Here, we only remind that 
A = 1 N eg) diag(l, 1, 1, 1, —I, —I, —I, —I) is the 8N x 8N 
supermatrix whose eigenvalues distinguish between the 
matrix blocks associated to the advanced and retarded 
Green functions and that J is a source matrix which is 
used to generate the Green function correlators by taking 
suitable derivatives with respect to its matrix elements 
II- 

Normally, H contains a certain number of gaussian 
random variables and it is possible to perform the en- 
semble average analytically which provides after a some 
transformations and approximations the non-linear su- 
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permatrix model p-|5||7]]. For the case of Cauchy band 
matrices, this approach has to be modified because of 
the non-gaussian distribution of its matrix elements. Our 
strategy will be to perform analytically by a new method 
the ensemble average with respect to the Cauchy dis- 
tributed diagonal matrix elements rjj while the coupling 
matrix U is kept fixed. This will provide a a model for- 
mulation for an arbitrary coupling matrix U . The en- 
semble average with respect to U can then be performed 
at a later stage when possible and convenient. As a first 
application of this procedure, we will derive an analytic 
expression for the localization length of the Cauchy band 
matrix model. 

The main idea is to replace the diagonal matrix ele- 
ments rjj by random phases e' l2ipj via r)j = (W/ir) tan^ . 
These phases are uniformly distributed on the unit cercle 
and can be formally identified with the scattering matrix 
of a chaotic cavity ideally coupled to one scattering chan- 
nel. Describing the Hamiltonian of the chaotic cavity by 
an M x M random matrix hj drawn from the gaussian 
orthogonal ensemble (GOE) one can apply the following 
substitution Jlyl : 



^ = (W/tt) tanv»j = (W/tt) A] hj' A 



(9) 



where Aj is an M -dimensional vector with components 
(Aj)i = 5i t i and hj is a real symmetric gaussian random 
matrix with variance ((hj)f k ) = (1 + Sik)/M. The substi- 
tution (||) is valid in the limit M — > oo which will be taken 
at the end of the calculations. We mention that this ap- 
proach is related to recent work concerning the derivation 
of the a model by performing phase averages in models 
involving unitary quantum maps for classically chaotic 
systems such as the kicked rotator or rough billiards 
pp| . The replacement (^) and the following steps follow 
very closely the method described in ref . Q . 

Inserting (|^) in the generating functional (||) and ap- 
plying suitable transformations inside the superdetermi- 
nant, one obtains an effective model defined on a Hilbert 
space of dimension N ■ M and characterized by an ef- 
fective Hamiltonian containing the GOE-matrices hj in 
its diagonal blocks and (—W/tt) Aj(D~ 1 )jk At in its off- 
diagonal jfc-block. Here D is the supermatrix appearing 
inside the superdeterminant of (j^) with rjj put to zero. 
This model is quite similar to the IWZ-model used in 
ref. ||] and following the standard procedure to derive 
the a model, we obtain (F(J)) V = J DQ e - c ^ with 
the action 



£(Q) = -Str ln(£ + (f + fe)A - U + i% ■ 



J 



(10) 



where Q contains in its diagonal blocks 8x8 superma- 
trices Qj belonging to Efetov's coset space for the or- 
thogonal symmetry class Q. We mention some impor- 
tant points concerning this result : (i) it is exact since 
the saddle-point approximation to derive the a model 



becomes exact in the limit M — > oo; (ii) its range of 
applicability goes far beyond the particular model con- 
sidered here since it is valid for an arbitrary coupling 
matrix U and its derivation only needs the Cauchy dis- 
tribution for the diagonal matrix elements; (iii) it gener- 
alizes the above mentioned replacement rule for one-point 
functions to the case of two-point functions : rjk + iOA — > 
i (W/tt) Qk where Qk is now a dynamical variable over 
which we have to integrate. 

To analyze the localization properties, we may put 
oj = and expand the action (10) for the long wave 
length limit. In the regime Uq b <C W, this gives apart 
from the source term the action 



C(Q) 
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dxUU) Str[(d x Q(x)f 



(11) 



with the quantity ^(U) = (4tt 2 /W 2 ) P E£ B+J . H 

£z(E/) were independent of x the action ( |ll| ) would corre- 
spond to the standard a model for quasi one-dimensional 
wires with the localization length £ = £ X (U) j|,[3]]. This 
is actually the case for gaussian distributed U xy for 
which £ X (U) w (£ x (U))u = 4ir 2 b 3 U§/3W 2 is a selfav- 
eraging quantity. However, for the Cauchy band ma- 
trix model £ X {U) is strongly fluctuating and its aver- 
age does not even exist. It is therefore necessary to 
take the s-dependence properly into account. A simi- 
lar situation has recently been encountered by Rupp et 
al. who derived the Id a model with an s-dependent 
"diffusion constant" from a hierachical random-matrix 
model for many-body states According to this 

work all quantities that can be extracted from the Id a 
model only depend on the rescaled length variable s(x) = 
Jo dy^y 1 ^). This justifies the intuitive average proce- 

dure T 1 = lim^i J Q L dx ^\U) « (^(U)). For- 
tunately, this average is finite and well defined for the 
Cauchy band matrix case (if b > 3). It can be easily 
evaluated by the same techniques used above for the lo- 
cal density of states [see eqs. (j§)-(0)]- For b ^> 1, we 
obtain : 



2tt b A UZ , 



7r t a\ 
16 VAV 



(12) 



where the dimensionless parameter 7 counts the number 
of well coupled levels inside a strip of typical length b, 
A = W/2b is the effective level spacing of such a strip, 
and a is the characteristic energy scale for the local den- 
sity of states (^|). (For the more general case of Levy 
distributed U jk , we find: £ - b 2+2 ^ 1 U$/W 2 .) The rela- 
tion (|l^) compares to £ ~ b (r/A) for the case of ordinary 
band matrices |Q,|l2|. However, it contradicts the expres- 
sion £ ~ &a/7 which was with proper translation of nota- 
tions numerically obtained from a similar Cauchy band 
matrix model in ref. We attribute this to the ex- 

treme numerical difficulty to access the regime 1 <ti 7 b 
where we expect (O) to be valid. For 7 < 1 we enter the 
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perturbative regime with £ ~ 6/111(7 1 ) while for 7 > b 
the localization length is likely to saturate at £ ~ b 2 . 




FIG. 1. The ratio f/6 as a function of the effective param- 
eter 7 ff. defined in eq. (]f^). The data points are numerical 
values, the dashed line corresponds to the limit £/& w 7 cff . 
for 7 c ff. S> f and the full line is the scaling curve (|f3|). The 
inset shows £/b versus 7. The different curves correspond to 
the values b = 10, 20, 30, 40, 50, 70, 85, 100 with b = 1 for the 
lowest curve. 

To verify this picture, we have therefore numerically 
determined the localization length by the recursive Green 
function method for 0.001 < 7 < 300 and 10 < b < 100. 
The inset of Fig. [j] shows the ratio £/& as a function of 
7. At first sight, it is indeed very difficult to verify the 
analytical result ( |l2| ) and a direct numerical fit (e. g. for 
b = 100 and 2 < 7 < 30) gives the power law £/b ~ 7 6 
rather than £/& w 7. To understand this, we note that 
the number of well mixed levels inside one bandwidth 
obviously saturates at b for 7 > b. It appears there- 
fore reasonable to replace in (|l^) 7 by an effective value 
7efi. that interpolates between 7 and b. Actually, we find 
that our numerical data can be well approximated by the 
scaling curve (see Fig. 

i „ 1-5 

b ~ ln(l + 1.5/ 7e ff.) 



7eff. 
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=- (13) 

[1 + v /3^t7W 

which reproduces the perturbative limit for 7 < 1, the 
analytical result ( |l2| ) for 1 -C 7 <C b and the behavior 
£ ~ b 2 for 7^6. Note that eq. ( |l3| ) provides for the 
case 7/6 <C 1 relative large corrections of order 1/7/6 ex- 
plaining the numerical difficulties to clearly identify this 
regime. 

In summary, we have obtained analytical results for 
the local density of states and the localization length for 
Cauchy random band matrices. We find that the Breit- 
Wigner form is replaced by a Levy distribution with in- 
dex /j,= l/2 and that the Breit-Wigner width T for ordi- 
nary band matrices is replaced by a new enhanced energy 
scale a ~ bT. This energy scale also determines the lo- 
calization length £ ~ 6 (a/A). From the technical point 
of view, we have derived a a model formulation for arbi- 
trary matrix ensembles with Cauchy distributed diagonal 
disorder. 

The author thanks G. Caldara, B. Georgeot, J. Lages 
and D. Shepelyansky for useful discussions. 
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